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Abstract 

One-dimensional  perturbations  of  a  two-component  neutral 
fluid  about  absolute  Maxwellian  equilibrium  states  (with  equal 
temperatures)  are  examined.   The  starting  equations  are  two 
Boltzmann  equations  coupled  through  diff erence-in-velocity 
terms.   A  Fourier-Laplace  transform  analysis  indicates  that 
the  resulting  modes  are  all  strongly  damped.   It  is  concluded 
that  a  momentum  exchange  mechanism  alone  cannot  support  the 
collective  sound  mode.   Asymptotic  forms  for  decay  frequency 
are  found  for  small  and  large  k  waves.   A  thermodynamic 
estimate  is  given  for  the  sound  speed  of  a  multi-component 
gas  in  which  the  number  of  degrees  of  freedom  of  each  component 
is  distinct. 
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Wave  Propagation  In  a  Multi-Component  Neutral  Fluid: 
The  Influence  of  Momentum  Exchange 

I.    Introduction  and  Summary  of  Results. 

When  investigating  the  fluctuations  associated  with  a 
multi-component  fluid  a  very  popular  course  of  action  used  to 
couple  the  components  is  to  use  momentum  exchange  terms  in  the 
equations  of  motion.   This  is  true  both  in  the  microscopic  vein 
[1,2,3],  where  the  formidable  collision  integrals  are  replaced 
with  simpler  diff erence-in-velocity  flow  terms,  and  it  is 
especially  true  in  the  macroscopic  approach  [4,5]  where  the 
individual  momentum  equations  are  supplanted  with  inhomogeneous 
momentum  exchange  terms.   This  latter  approach  has  recently  been 
extended  to  analysis  of  partially  ionized  gases  [6-11]  (i.e., 
three-component  fluids)  with  specific  regard  to  problems 
pertaining  to  ionospheric  propagations. 

In  all  of  these  problems  an  important  goal  is  that  of 
recapturing  the  multi-component  sound  speed.   It  is  the  purpose 
of  this  paper  to  show  that  any  analysis  of  wave  propagation  in 
a  multi-component  fluid  whose  only  interaction  mechanism  is  one 
of  momentum  exchange  will  exhibit  only  strongly  damped  modes, 
and  consequently  can  not  support  the  proper  (undamped  in  one 
dimension)  multi-component  sound  phenomenon.  This  observation, 
of  course,  renders  the  included  model  unrealistic. 
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The  Included  starting  equations  are  two  one-dlmenslonal 
Boltzmann  equations,  which  have  been  linearized  about  absolute 
Maxwelllan  state  with  equal  temperatures.   These  equations  are 
coupled  through  dlff erence-ln  -velocity  forcing  terms.   A 
Fourier-Laplace  transform  analysis  Is  applied  to  a  generalized 
boundary-Initial  value  problem.   The  results  are  obtained  by 
examining  the  poles  of  the  Integrand  of  the  Inverse  transform. 
These,  In  turn.  Indicate,  that  there  are  no  growing  modes 
(i.e.,  the  only  solutions  are  for  Re  s  <  0,   where   s   Is  the 
time  transform  variable).   Analytic  continuation  of  the  transform 
also  yields  no  roots,  both  for  Re  s  =  0   (pure  propagation)  and 
In  a  neighborhood  Inflnlteslmally  close  to  the  Imaginary   s 
axis  (propagation  with  small  attenuation).   The  only  solutions 
that  remain  are  those  which  are  strongly  damped.   Two  specific 
extremes  are  calculated.   The  first  pertains  to  large  k 

(k  =  wave  number)  phenomena.   These  waves  decay  with  a  frequency 

2 
which  varies  as  k  .   The  second  example  pertains  to  small  k 

phenomena,  whose  decay  frequency  varies  In  a  slow  logarithmic - 

like  manner  with  k. 

In  retrospect.  It  Is  clear  that  a  momentum  exchange 

mechanism  will  not  support  the  sound  phenomena.   This  mechanism 


* 

This,  of  course.  Is  equivalent  to  the  fact  that  there  are 
no  "normal"  modes. 
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serves  primarily  to  quench  macroscopic  velocity  differences  but 
cannot  force  local  equilibrium.   Such  equilibrium  is  attained 
In  a  one-component  fluid,  for  Instance,  through  a  BGK-llke  [12] 
interaction. 

As  pointed  out  above,  in  studies  of  wave  propagation  in 
partially  ionized  gases,  the  mechanism  of  momentum  exchange  is 
used  to  couple  the  three  components.   However,  Inasmuch  as  two 
of  the  species  are  charged,  there  is  another  coupling  mechanism, 
namely  that  afforded  by  Maxwell's  equations.   This  electrical 
coupling  enforces  an  equilibrium  corresponding  to  the  state  of 
charge  neutrality.   So  it  is,  for  instance,  that  Landau's  [lj>] 
analysis  which  Includes  only  long  range  averaged  electric  fields, 
without  any  short  range  collisions,  gives  propagating  waves, 
which  in  the  limit  of  the  long  wavelengths  severely  resemble 
sound  phenomena.   However,  the  true  fluid  dynamical  long  wave- 
length sound  mode  (e.g.,  as  appears  in  the  Lundqulst  equations) 
has  been  shown  to  vanish  in  the  limit  of  vanishing  short  range 
collision  frequency  [^] . 

A  recent  macroscopic  analysis  by  Tanenbaum  and  Mlntzer  [ll] 
of  wave  propagation  in  a  partially  ionized  gas  includes  both 
frlctional  and  electrical  coupling  between  the  species.   In 
addition,  equations  of  state  are  Included  for  each  of  the  three 
components,  which  in  turn  Insures  Internal  equilibrium  for  the 
separate  species.   Under  these  conditions  the  momentum  exchange 
mechanism  may  well  enforce  total  equilibrium  Inasmuch  as  it  Is 
driving  together  components  already  in  equilibrium.   Indeed,  in 


the  limit  of  low  frequencies,  the  adlabatlc  sound  speed,  \/yp/p ' 
is  recaptured,  where   p  and  p   are  equilibrium  total  pressure 
and  mass  density  respectively. 

A  brief  thermodynamic  derivation  of  the  cooperative  sound 
speed  Is  Included  In  the  appendix. 


II.   Analysis 

Our  starting  equations  are  two  one-dlmenslonal  augmented 
Boltzmann  equations  [3]  for  the  perturbation  distribution   f: 


(1)        if-5i= voij«^- 


(2)         i-?if=-vot«'^" 


where 

( 5  )  Au  =  u  -  u 


(^) 


mn  mn 

[X     —      AA 


mn^  +  mn 
c-     o 


(5)  KT  =  m  C^  -  m  C^ 


(6)  f^  =  \-   exp(-eV2c2), 

°    (27^)2  C 
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Furthermore,  the  distributions  are  normalized  so  that. 


(7)  /f  =  n 

V  I  /  •'00 


(8)  /fl  -  n^u. 


The  variables  associated  with  "gas  No.  -l"  are  roofed'  while 
those  of  "gas  No.  2"  are  bare.   The  two  gases  have  the  same 
temperature   T   in  the  equilibrium  state  which,  in  turn,  is 
specified  by  the  absolute  Maxwellians   f   and  f  .   The  perturbed 
state  away  from  equilibrium  is  given  by  the  total  distribution. 


(9)  P  =  f  +  f :     P  =  f  +  f 

^-^ '  o  o 


Number  density  is   n,   macroscopic  velocity  is   u,   m  is  mass, 

M,  is  reduced  mass,   C   is  thermal  speed,   ^   is  microscopic  speed, 

and   V   is  cross-collision  frequency. 

The  problem  that  is  considered  is  a  mixed  initial-boundary 
value  problem  with   f(x,0)  =  h(x)   and  f(t  00,0)  =  0.   The 
method  of  solution  is  that  of  the  Laplace-Pourier  transform: 


00   00 
r 


(10)     f*(s,k)  - 


dx 


,,   -st+ikx  „/,  ^\ 
dt  e        f(t,xj 


-00    o 


(Re  s  >  0) 


■X-  ^ 

All  distributions  are  functions  of  velocity  %. 

More  accurately^  M.  is  a  density  weighted  reduced  mass, 
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with  corresponding  inverse. 


(11)  f(t,x)   = 


i(27r)' 


00  7+ic 

r       ^ 


dk 


ds   e^^-^^^  f*(s,k) 


'-00        y-ioo 


Multiplying  equations  (l,2)  through  by  the  Integral  operator 
/  dx  /  dt  exp[-st+ikx]   gives  the  desired  transformed  equation. 


(12) 


M-  e  Au 


[s  -  ikllf  -  h(x)  =  v^  ^  ^  ^  f^ 


(13) 


[s  -  ikllT  -  h(x)  =  -v^  ^  ^  ^  f^ 

m   c 


* 


Solving  for   f   gives 


(14) 


^      (n/m)  e  (AuVc^)f 
f   =  v^- 


s-ik| 


o  ^  h(x 


s-ik£  ' 


with  a  similar  expression  for  f  .   Forming  the  velocity  moment 
gives 


(15) 


*    ^]_M-   ^^* 


u   =  - 


imkC 


[u  -u] 


(27r)"^e-^  /^  v^dv 
z-v 


+  <h> 


—  00 


-  9  - 


00 


(16) 


u   = 


1     ^^ 


Imlcc 


+  <h> 


(27r)  ^e  '  ^  ^  v^dv 


z-v 


J 

—  oo 


z  =  s/lCk 

where  <h>   Is  the  Implied  Initial  distribution  moment.   Sub- 
traction of  these  two  equations  gives  one  equation  for  the 
single  unknown  Au  , 


(17) 


Au 


A  <h> 


^  (Cl(z)  +  Cl(z))  +  1 


where 


(18) 


a  =  KT^v^n 


(19) 


r 


I(z)  = 


(2Tr)' 


2  -V  /2  -, 
V  e   "^   dv 

z-v 


-00 


(20) 


r 


<Ah>  =  C 


2  /o  ^ 

e"^  /^[h-h]vdv 

z-v 


kJ 


The  inverse  of  (l?)  appears  as 
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(21) 


Au 


i2iry 


00 

r 


dke 


\j 


Y+ioo 


•Ikx 


dse^^A<h> 


[Cl(z)  +  Cl(z)] 


lak 


+  1 


y-lc 


In  the  event  of  vanishing  coupling  a  — ^  oo,   and  only  the 
A<h>   factor  remains  which.  In  turn,  Induces  the  streaming  solution 
Au  =  Au(0).   More  generally,  both  the  poles  of   A  <h>   and  the 
roots  of 


(22) 


/\     /\ 


Cl(z)  +  Cl(z 


■lak 


contribute  to  the  time  development  of  Au.   In  uncovering  the 
roots  of  eq.(22)  It  Is  Important  to  recall  that  the  definition 
of  f  ,  and  therefore  Au  ,   Includes  Re  s  >  0,   which  (for 
positive  k)   relates  to  Im  z  <  0. 

Examination  of  the   I  functions  quickly  Indicates  that 
there  are  no  roots  for   Im  z  <  0  which,  In  turn,  relates  to  the 
fact  that  the  starting  equations  (l,2)  with  the  Included  boundary- 
Initial  conditions  do  not  Induce  growing  modes.   A  separate 
analysis  for   z   on  the  real  axis  Indicates  that  there  are  also 
no  roots  there  —so  that  there  are  no  purely  propagating  modes. 
Indeed,  continuing  the   I  functions  Into  the  upper  half  z-plane 

gives  no  roots  In  an  open  region  above  the  real  axis so  that 

there  are  no  propagating  modes  which  are  only  slightly  damped. 
The  only  solutions  that  equation  (2l)  has  are  In  the  remaining 
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section  of  the  upper  half  z-plane. 

In  order  to  exhibit  these  facts  we  first  rewrite 
l(z)  In  the  form. 


(19  )     I(z) 


(27r)' 


V 


e-  /^  vg(z-v)dv     (i„,,o), 

lz-v| 


so  that 


(23)       Im  I(z)  =  D  Im  z  =  -Dy  <  0 


(y  >  0), 


Where   D  Is  the  Implied  real  positive  definite  function. 
For   z   approaching  the  real  axis  the  contour  must 
be  distorted  as  depicted  in  Pig.  1. 


Z-plane 


X-€ 


X+€ 


Figure  I.  Continuation   onto  the   Real    Axis. 
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In  the  limit  as   e  7  0,   one  obtains  the  proper 

analytic  combination  of  l(z)   Into  the  real  axis  [13]^  which 
appear  as. 


r 


(24)      I(x) 


(27r)2 


Vj 


s:l^i^!av,i(,/2)i,v-'/2. 


x-v 


(Im  z  =  0) 


where  the   P:   symbol  Indicates  that  the  principal  part  Is 
to  be  taken.   It  follows  that  along  the  real  axis. 


(25) 


Im  I(x)  =  (Tr/2)^  x^  e 


\      2   -x2/2 


a  positive  definite  form. 

In  the  upper  half  plane  the  contour  Is  distorted  as 
depicted  In  Figure  2, 


0 


z-plane 
Im  z  =  0 


Figure  2.        Continuation  into  the  Upper  Half  Plane. 
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to  obtain 


(26) 


I(z) 


(27r)- 


r 


-00 


Z-V 


+  l(7r/2)2  z^e"^  /2 
(Im  z  >  0), 


or  equlvalently. 


(27) 


I(x+iy) 


(27r)' 


'J 

(¥-v)v^  e-'^  /^  dv 
(x-v)  +  iy 


—  00 


7^^i  /„2  _.2,^...__^  .  -  [x^-y^+2ixy  ]  ^ 


+  1(|)2  (x^-y^+21xy)  e 


and  also 


(28) 


I(z)  = 


(27r)' 


/-   ^  2   -V  /2 
(z-vjv   e   ' 

2 


z-v 


dv  ,  .  /ttnI   2  -z  /2 
—  +  ligj^  z  e   ^  . 


>J 


a.   Growing  Modes 

We  are  now  prepared  to  exhibit  formally  the  statements 
made  above.   First,  for   Im  z  <  0,   comparison  of  equations 
(21-23)  Immediately  reveals  that  there  can  be  no  solutions  In 
this  domain,  granted  that  k  >  0.   This,  of  course.  Implies  that 
there  are  similarly  no  growing  modes  for  k  <  0  ,  inasmuch  as  the 
starting  equations  are  invariant  under  reflexion  of  coordinates. 


_  1^  - 


I 

However,  it  Is  also  evident  If  we  merely  set   k  ^  k  =  -k, 

so  that   z  — ^  z  =  -z,   and   Im  I(z)   tiecomes  positive,  and 
equation  (22)  again  has  no  solutions  in  the  domain  of  growing 
modes. 


b.  Purely  Propagating  Modes 

The  purely  propagating  waves  correspond  to  roots  along  the 
real  z-axis.   Comparison  of  the  dispersion  equation  (2l)  and 
equation  (24)  again  shows  that  there  are  no  solutions  for  k  >  0. 
The  above  space  reversal  argument  shows  that  there  are  also  no 

solutions  for  k  <  0.   More  explicitly,  consider  again  that 

I  1  ,  , 

k  ^  k  =  -k,   so  that   z  ^  z   =  -z.   Inasmuch  as   I(z)   is 

defined  primarily  for  Re  s  >  0,   the  real  axis  must  now  be 

approached  from  above  and  the  Cauchy  contribution  reverses  sign 

since  the  direction  of  the  loop  integral  is  reversed.   It 

follows  that  the  imaginary  parts  of  both  the  right  and  left 

sides  of  (21)  change  sign,  and  solution  is  still  impossible  along 

the  real  axis. 

c.  Slightly  Damped  Modes 

These  modes  are  those  lying  in  the  region  of  the  z-plane 
depicted  in  Figure  3.   The  proof  that  there  are  no  slightly 
damped  modes  is  given  first  for  the   e-wedged  region,  and  then 
for  the  region  about  the  origin. 


-15- 


Figure  3. 
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The  wedged   region   Is   given  by     0  <  0   <   e      and     tt   -   e  <  0   <  it. 
In   that   region,    sin  9   =  0(e)      and      cos  6   =   0(l). 
Substituting      z  =  re^^      Into  eq.    (28)   and   retaining  terms 
to   order  unity      gives 


(29) 


Im   I(z)    =   B(cos    2e)e"''   ^°^   ^    +   0(e), 


More  precisely,  the  strongly  damped  modes  lie  in  a 
cone  about  the  imaginary  z-axls  with  apex  Tr/4.   The 
estimate  given  in  Eq.(29)  begins  to  fall  drastically 
near   e  =  7r/4,  where   cos  £  ^  sin  e. 
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(where   B  Is  real  and  positive)  which,  in  turn,  indicates  that 
the  imaginary  part  of  the  dispersion  relation  equation  (22) 
cannot  he  satisfied  in  this  domain. 

In  the  remaining  domain  about  the  origin  (excluding  the 
imaginary  axis),  if  z  =  e  is  substituted  into  eq.  (26)  then 
one  obtains. 


(30)      I(z)  =  -  — ^ 

(27r)2 


V  e-^^/2^v[l  +  (|)  +  (f)2  +  ...] 


—  00 


+  i(|)^  z2(l  -  ^  +  ...), 


so  that  to  terms  of   0(z  ), 


'Z^  .2  ^  ^/  3 


(31)  I(z)  ~  -  z  +  i(|)  z^  +  0(z-^). 


Consider  now  the  real  part  of  the  left  side  of  the 
dispersion  equation  (2l).   To  terms  of   0(z),  there  results 


/>>     /N 


(32)  Re[Cl(z)  +  Cl(z)]  =  -  [Cx  +  Cx], 

which  cannot  be  zero  for  x  >  0,   so  that  eq.  (22)  has  no 
solutions  in  this  domain. 
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d.  Remaining  Solutions 

The  above  analysis  Indicates  that  the  only  solutions 
of  the  dispersion  equation  (22)  lie  in  the  sector  of  the 
upper  half  z -plane  shown  in  Fig.  4. 


Figure  4. 


Domain  of  Solution 


That  there  are  any  solutions  at  all  follows  from  the  entire 
nature  of  the   I   functions  about  the  point  at  infinity. 
That   I   is  entire  may  be  seen  by  relating  it  to  another 
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entire  function,  viz.,  the  Landau  [13]  function  L, 
(33) 


00 


L 


.  1 

-• T- 

(27r)2 


U 


-v^/2  ^ 
V  e   '^   dv 


v-z 


This  is  related  to  the   I   integral  through 


(34) 


-I  = 


(27r)' 


V  /2 
v(v-z+z)e   "^ 


dv 


v-z 


zL 


-00 


so  that  I  is  entire,  granted  that  L  is  entire  . 

We  will  now  exhibit  two  explicit  solutions  in  this 
remaining  domain,  both  of  which  lie  along  the  imaginary  axis, 

First  for   z  =  iy,   y  «  1,  [k  »  s/C] 
the  form  (26)  gives 


(35) 


i(iy)  = 


(27r)' 


-v^/2   2, 
e   '^   V  dv 


^(1L\\     2  ^yV2 


ly-v 


-  l(^) 


y  e' 


n- 


AiL 


(27r)^ 


^-v72   2 
e   ^   V 

2,  2 
y  +v 


Kf)*  y^e^'/s 


ly+o(y  ). 


\^ 


*  Indeed,  a  disperson  equation  identical  to  (22),  only  in  terms 
of  L  rather  than   I,  results  if  the  difference  between  the 
two  continuity  equations  are  combined  with  the  number  density 
moments  of  f  (instead  of  velocity  moments)  to  again  obtain 
one  equation  for  Au. 
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If  this.  In  turn.  Is  substituted  Into  (21)  there  results. 


(36)  -  l[Cy  +  Cy]  =  +if^=  -  — -° 


Where  we  have  recalled  that   z  =  ly  =  s/ick,   so  that   y  =  -s/Ck. 
The  root  of  (36)  appears  as 


(37)  s  =  -  k^   ° 


2v^M. 


Our  second  example  also  involves  roots  along  the  imaginary 
axis   z  =  iy;  however,  in  this  case   y  »  1.   The  Cauchy 
contribution  dominates  over  the  integral  part  in  this  region  so 
that  (26)  yields 

(38)  l(z)  ~  -  1(f)*  y2  ey^/2. 
The  dispersion  equation  (22)  now  reads, 

(39)  0(y)  =  (|)i  Ty^c  e^'/^  +  ?C   e^'/^  |  =  ka  =  -  |  (§). 

Since   y  =  y[c/c],  the  left  side  of  this  equation  may  be  written 

purely  as  a  function  of  y,  viz.,   G(y).   This  function  is  more 

rapidly  increasing  than  exponental  (for  large   y).   However,  to 

gain  a  rough  estimate  of  the  behavior  of  s  =  s(k),   let  us 

2  11 

suppose  that   G(y)  ~  e'^  ,  so  that  y  ~  -gn(ka)2  and  s  ~  ka  in(ka)2 
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The  actual   G  is  a  more  rapidly  increasing  function  of  y   so 
that,  in  truth,   s   varies  even  more  slowly  than  this  logarithmic 
estimate  indicates. 
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Appendix. 

The  formula  for  the  sound  speed  of  a  multl -component  gas, 
which  follows  from  fluid  dynamics  [l^].  Is 


2  at 


(A.l)  a^=|£. 


evaluated  at  constant  entropy  and  at  equilibrium.   The  total 
pressure   p   and  total  mass  density  p   are  related  to  the 
component  values   p.   and   p.,   respectively,  through 


(A. 2)  P  =  2  p. 


(A. 3)  p  =  2  p^. 

A  similar  relation  holds  between  the  total  Internal  energy 
E  and  the  component  energies   E. ,  viz.. 


(A. 4)  E  =  2  E. , 


The  equilibrium  state  of  each  component  obeys  the  relations. 


f^      r.      ^     .r  f^ 


(A. 5)      E.  =  ^  p.R.T.V  =  j—   p.V, 
(A. 6)      p.  =  p.R.T. , 


(A. 7)      T.  =  T   . 
'       1     o 
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The  mass  density  p.   Is  related  to  the  total  volume 

through  the  elemental  mass   m.   and  the  number  of  particles  in 

th 
the   i     species   N.   according  to 

m.N.    c. 
(A. 8)  p.  =_^  =  ^, 


(A. 9)  p  =    Zp^   =  ^  Zc^   =  ^    , 


where   c   Is  the  Implied  constant.   The  number  of  degrees  of 

th  1 

freedom  of  an  element  of  the   1  '   species  Is   f  .   The  constant 

R   Is   K/m.   where   K  Is  Boltzmann ' s  constant. 

In  order  to  calculate   Sp/Sp   at  constant  entropy  S,   we 

employ  the  first  law  of  thermodynamics,  which  under  the  said 

constraint  appears  as, 

(A. 10)         Tds  -  0  =  dE  +  pdV. 

Substituting  the  above  equilibrium  relations  gives. 


<-'  r-(fe-^^;)(^^ 


V 


If  we  assume  that  the  term  in  brackets  is  slowly  varying  in 
the  adiabatlc  variation  (A.IO),  then  (A.ll)  is  readily  integrated 
to  obtain  the  familiar  form, 

(A. 12)  p  =  ApT', 
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where  we  have  set  y     equal  to  the  bracketed  term  In  (A.ll),  and 
eliminated  V  through  equation  (A.9).   The  constant   A   is  as 
Implied.   Equation  (A. 12)  yields  the  desired  form: 


(A. 13)  a^  =  t(p/p). 


If  all  the   f   are  equal  then  the  sound  speed  becomes 

(A.14)       [^^]  ^=  {^^)  kt  (^ ,, 

^     '  f    2p .    ^  f  '    o    2n .  m . 

where   n   is  the  number  density  of  the   i    species, 
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